3D+1 Lorentz type soliton in air 



Lubomir M. Kovachev 
Institute of Electronics, Bulgarian Academy of Sciences, 
Tzarigradcko shossee 72,1784 Sofia, Bulgaria, e-mail:lubomirkovach@y ahoo.com 

(Dated: January 16, 2013) 

Up to now the long range filaments have been considered as a balance between Kerr focusing and 
defocusing by plasma generation in the nonlinear focus. However, it is difficult to apply the above 
explanation of filamentation in far-field zone. There are basically two main characteristics which 
remain the same at these distances - the super broad spectrum and the width of the core, while the 
power in a stable filament drops to the critical value for self-focusing. At such power the plasma and 
higher-order Kerr terms are too small to prevent self-focusing. We suggest here a new mechanism 
for stable soliton pulse propagation in far-away zone, where the power of the laser pulse is slightly 
£SJ ' above the critical one, and the pulse comprises super-broad spectra. For such pulses the diffraction 

is not paraxial and an initially symmetric Gaussian pulse takes parabolic form at several diffraction 
lengths . The stable soliton propagation appears as a balance between the divergent parabolic type 
diffraction of broadband optical pulses and the convergent nonlinear refractive index due to the 
intensity profile. We investigate more precisely the nonlinear third order polarization, using into 
account the carrier-to envelope phase. This additional phase transforms the third harmonic term 
. to THz or GHz one, depending on the spectral width of the pulse. 

PACS numbers: 42,65.-k, 42.65.Tg 
Keywords: analytical three dimensional bright soliton, supercontinuum generation 



OS 



I. INTRODUCTION 
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In the process of investigating the filamentation of a power femtosecond (fs) laser pulse many new physical effects 
have been observed, such as long-range self-channeling H - 3[, coherent and incoherent radial and forward THz emission 
[3-[H , asymmetric pulse shaping, super-broad spectra [7 - TJ H3| and others. The role of the different mechanisms in 
i-C ' near zone (up to 1 — 2 m from the source) has been investigated experimentally and by numerical simulations, and 
most processes in this zone are well explained (l2H15j . When a fs pulse with power of several P cr = 7r(0.61Ao) 2 / (8nrjn2) 
starts from the laser source a slice- by-slice self- focussing process takes place 16]. At a distance of one-two meters the 
pulse self-compresses, enlarging the k z spectrum to super-broad asymmetric spectrum Ak z s» k . The process increases 
^ ■ the core intensity up to tens of 10 13 W/cm 2 , where different types of plasma ionization, multi-photon processes and 
higher-order Kerr terms appear [l7^ . Usually, the basic model of propagation in near the zone is a scalar spatio- 
temporal paraxial equation including all the above mentioned mechanisms [T2I [ill . [l7j . The basic model is natural in 
the near zone because of the fact that the initial fs pulse contains a narrow-band spectrum Ak z << kg. Thus, the 
paraxial spatio-temporal model gives a good explanation of nonlinear phenomena such as conical emission, X-waves, 
^Tj ■ spectral broadening to the high frequency region and others. In far-away zone (propagation distance more than 2 — 3 
meters) plasma ionization and higher-order Kerr terms are admitted also as necessary for a balance between the 
I self-focussing and plasma defocussing and for obtaining long range self-channeling in gases. 
• • ' However, the above explanation of filamentation is difficult to apply in far-away zone. There are basically two main 
characteristics which remain the same at these distances - the superbroad spectrum and the width of the core, while 
k> the intensity in a stable filament drops to a value of 10 12 W/cm 2 [H, U?]]. The plasma and higher-order Kerr terms 
are too small to prevent self-focussing. The observation of long-range self-channeling [l8^20] without ionization also 
_ leads to change the role of plasma in the laser filamentation. 

In addition, there are difficulties with the physical interpretation of the THz radiation as a result of plasma 
generation. The plasma strings formed during filamentation should emit incoherent THz radiation in a direction 
orthogonal to the propagation axis. The nature of the THz emission, measured in [6j is different. Instead of being 
emitted radially, it is confined to a very narrow cone in the forward direction. The contribution from ionization 
in far-away zone is negligible fl7| and this is the reason to look for other physical mechanism which could cause 
THz or GHz radiations. Our analysis on the third order nonlinear polarization of pulses with broadband spectrum 
indicates that the nonlinear term in the corresponding envelope equation oscillates with frequency proportional to the 
group and phase velocity difference fL,; = 3(kov p h — v gr Ak z ). Actually, this is three times the well-known Carrier-to 
Envelope Phase (CEP) difference [21] . This oscillation induces THz generation, where the generated frequency is 
exactly VLthz = 93GHz for a pulse with superbroad spectrum Afc z f=a fco with carryier wavelength 800 nm. 

Physically, one dimensional Schrodinger solitons in fibers appear as a balance between the Kerr nonlinearity and 
the negative dispersion [22H241 ] . On the other hand, if we try to find 2D+1 and spatio-temporal solitons in Kerr 
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media, the numerical and the real experiments demonstrate that there is no balance between the plane wave paraxial 
diffraction - dispersion and the Kerr nonlinearity. This leads to instability and self-focusing of a laser beam or initially 
narrow band optical pulse. Recently Serkin in [25| suggested stable soliton propagation and reducing the 3D soliton 
problem to one dimensional, with introducing trapping potential in Bose - Einstein condensates. 

In this paper we present a new mathematical model, on the basis of the Amplitude Envelope (AE) equation, up to 
second order of dispersion, without using paraxial approximation. In the non-paraxial zone the diffraction of pulses 
with superbroad spectrum or pulses with a few cycles under the envelope is closer to wave type 26]. For such pulses, a 
new physical mechanism of balance between nonparaxial (wave-type diffraction) and third order nonlinearity appears. 
Exact analytical three-dimensional bright solitons in this regime are found. 



II. LINEAR REGIME OF NARROW BAND AND BROAD BAND OPTICAL PULSES 



The paraxial spatio-temporal envelope equation governs well the transverse diffraction and the dispersion of fs 
pulses up to 6 — 7 cycles under the envelope. This equation relies on one approximation obtained after neglecting the 
second derivative in the propagation direction and the second derivative in time from the wave equation 2jJ or from 
the 31? + 1 AE equation 28j. In air, the series of k 2 (w) are strongly convergent up to one cycle under the envelope 
and this is the reason why the AE equation is correct up to the single-cycle regime. 

The linearized AE, governing the propagation of laser pulses when the dispersion is limited to second order, is: 
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where f3 = k^kav^ is a number representing the influence of the second order dispersion. In vacuum and dispesionless 
media the following Diffraction Equation (DE) (v ~ c ) is obtained: 
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We solve AE (p} and DE ^ by applying spatial Fourier transformation to the amplitude functions A and V. The 
fundamental solutions of the Fourier images A and V in (k x ,k y ,Ak Zl t) space are: 
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A = A(k x , k y , Ak z ,t = 0) x 
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V = V(k x ,k y , Ak z , t = 0) exp <iv[k ± Jk x 2 + k y z + {Ak z - k Q ) 2 t 
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respectively. In air j3 ~ 2.1 x 10~ 5 , AE (JXJ) is equal to DE @, and the dispersion is negligible compared to the 
diffraction. We solve analytically the convolution problem ((4]) for initial Gaussian light bullet of the kind V(x, y,z,t — 
0) = exp (—(a; 2 + y 2 + z 2 )/2tq). The corresponding solution is: 
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where f — \f x 2 + y 2 + (z — ir^fco) 2 . On the other hand, multiplying the solution ([S]) with the carrier phase, we obtain 
solution of the wave equation E (x, y, z,t) — V (x, y, z, t) exp (i(k$z — uj^i)), where u>o and fco are the carrier frequency 
and carrier wave number in the wave packet: 
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A systematic study on the different kinds of exact solutions and methods for solving wave equation © was per- 
formed recently in (2||. Here, as in [26| we suggest another method: Starting with the ansatz E(x,y,z,t) = 
V (x, y, z, t) exp (i(koZ — wot)), we separate the main phase and reduce the wave equation to 3D + 1 parabolic type 
one ©. Thus, the initial value problem can be solved and exact (or numerical) solutions of the corresponding 
amplitude equation @ can be obtained. The solution ([5]), multiplied by the main phase, gives an exact solution ([7]) 
of the wave equation (|6]). To investigate the evolution of optical pulses at long distances, it is convenient to rewrite 
AE (JT|) equation in Galilean coordinate system t' = t; z' = z — v gr t: 
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Pulses governed by DE ^ move with phase velocity and the transformation is t' = t; z' = z — vt: 
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Here, = + j^p denotes the transverse Laplace operator. The corresponding fundamental solution of AE 
equation © in Galilean coordinates is: 



A G (k x ,k y ,Ak Zl t) = A G (k x ,k y ,Ak z ,t = 0) x 
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while the fundamental solution of DE (|9]) becomes: 



V G = V G (k x ,k y ,Ak z ,t = 0)x 
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The analytical solution of (fTTj) for initial pulse in the form of Gaussian bullet is the same as ([5]), but with new radial 
component f = x 2 + y 2 + (z + vt — ir^ko) 2 translated in space and time. The numerical and analytical solutions of 
AE (JIJ and DE ([2]) are equal to the solutions of the equations AE ([8]) and DE (J9j) in Galilean coordinates with only 
one difference: in Laboratory frame the solutions translate in z-direction , while in Galilean frame the solutions stay 
in the centrum of the coordinate system. 

The basic theoretical studies governed laser pulse propagation have been performed in so called "local time" 
coordinates z — z\ t = t — z/v gr . In order to compare our investigation with these results, we need to rewrite 
AE equation ([!]) for the amplitude function A in the same coordinate system. Thus Eq. ([1} becomes: 
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Since this is a parabolic type equation with low order derivative on z, we apply Fourier transform to the amplitude 
function in form: A (k x , k y , Aw, z) — FFF [A (x,y, z, t)], where FFF denotes 3D Fourier transform in x,y,r space 
and Aw = w — ujq; Ak z = Auj/v gr are the spectral widths in frequency and wave vector domains correspondingly. 
The following ordinary differential equation in (k x , k y , Aw, z) space is obtained: 
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As can be seen from (|13[) . if the second derivative on z is neglected, then the paraxial spatio-temporal approximation 
is valid. Equation (fT3"|) is more general and we will estimate where we can apply spatio-temporal paraxial optics (PO), 
and where PO does not works. The fundamental solution of (fT3")) is: 
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The analysis of the fundamental solution (fT4|) of the equation (1131) is performed in two basic cases: 

a: Narrow band pulses - from nanosecond up to 50 — 100 femtosecond laser pulses, where the conditions: 
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are satisfied, and the wave vector's difference ko — Au!/v gr can be replaced by fco. Using the low order of the Taylor 
expansion and the minus sign in front of the square root from the initial conditions, equation (1141) is transformed in 
a spatio - temporal paraxial generalization of the kind: 



A(k x ,k y ,Au,z) = A(k Xl k y , Aw, 0) exp 

From (1161) the evolution of the narrow band pulses becomes obvious: while the transverse projection of the pulses 
enlarges by the Fresnel's law, the longitudinal temporal shape will be enlarged in the same away, proportionally to the 
dispersion parameter /3. Such shaping of pulses with initially narrow band spectrum is demonstrated in Fig.l, where 
the typical Fresnel diffraction of the intensity profile (spot (x, y) projection) is presented. The numerical experiment 
is performed for 100 femtosecond Gaussian initial pulse at A = 800 nm, Ak z << fco, zq = 30/im, r^{x,y) — 60/im, 
with 37.5 cycles under envelope propagating in air (ft = 2.1 x 10~ 5 ). The result is obtained by solving numerically 
the inverse Fourier transform of the fundamental solution (ITU of the AE equation in the local time frame (Ti"2"j) . The 
spot enlarges twice at one diffraction length Zdiff — r 2 ko. Fig. 2 presents the intensity side (x,t) projection of the 
same pulse. We should note that while the spot ((x,y) projection) enlarges considerably due to the Fresnel law, the 
longitudinal time shape (the r projection) remains the same on several diffraction lengths from the small dispersion 
in air. The diffraction - dispersion picture, presented by the side (x, r) projection of the pulse, gives idea of what 
should happen in the nonlinear regime: the plane wave diffraction with a combination of parabolic type nonlinear 
Kerr focusing always leads to self- focusing for narrow-band (Ak z « fco) pulses. The same Taylor expansion for 
narrow band pulses can be performed to fundamental solutions of the equation in Laboratory §3§ and Galilean (|10[) 
frames. 

b: broad band pulses - from attosecond up to 20 — 30 femtosecond pulses, where the conditions: 
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z=0 z=U*$12 z=Zdiff z = 3z dlf /2 

FIG. 1: Plot of the waist (intensity's) projection |A(x,y)| 2 of a 100 fs Gaussian pulse at A = 800 nm, with initial spot 
ro — 60 /im, and longitudinal spatial pulse duration zo = 30 am, as solution of the linear equation in local time (|12[) on 
distances expressed by diffraction lengths. The spot deformation satisfies the Fresnel diffraction law and on one diffraction 
length z = Zdiff the diameter of the spot increases twice, while the maximum of the pulse decreases with the same factor. 




z=0 z=ZdiFf''2 z-zas z-3zdiff/2 



FIG. 2: Side (x,t) projection of the intensity |A(:e,t)| 2 for the same optical pulse as in Fig. 1. The (x,y) projection of the 
pulse diffracts considerably following the Fresnel law, while the (r) projection on several diffraction lengths preserves its initial 
shape due to the small dispersion. The diffraction - dispersion picture, presented by the side (a;, r) projection, gives idea of 
what should happen in the nonlinear regime: the plane wave diffraction with a combination of parabolic type nonlinear Kerr 
focusing always leads to self-focusing for narrow-band (Ak z « fco) pulses. 

are satisfied. In this case we can not use Taylor expansion of the spectral kernels in Laboratory ([3]), Galilean (fTU)) 
and local time ((14]) frames. The spectral kernels are in square root and we can expect evolution governed by wave 
diffraction. That why for broadband pulses we can expect curvature (parabolic deformation) of the intensity profile of 
the (x, z) or (x, r) side projection. Fig 3. present the evolution of the intensity (side (x, z) projection) of a normalized 
10 fs Gaussian initial pulse at A = 800 nm; Ak z ~ fco/3; zq = ro/2; and only 3 cycles under the envelope (broadband 
pulse), obtained numerically from AE equation (jHJ in Galilean frame. The solution confirms the experimentally 
observed parabolic type diffraction for few cycle pulses. And here appears the main physical question for stable pulse 
propagation in nonlinear regime: Is it possible for the divergent parabolic intensity distribution due to non-paraxial 
diffraction to be compensated by the converged parabolic type nonlinear Kerr focusing? If this is the case, then a 
stable soliton pulse propagation exists. As we show below, only for broadband pulses one-directional soliton solution 
of the corresponding nonlinear equations can be found. 




z-L) z=z<ijfi z=2zdifl z-3zdiff 



FIG. 3: Side (x, z) projection of the intensity \A(x, z')\ 2 for a normalized 10 fs Gaussian initial pulse at A = 800 nm, Ak z ~ fco/3, 
zo = ro/2, and only 3 cycles under the envelope (large-band pulse Ak z « fco), obtained numerically from the AE equation 
in Galilean frame. At 3 diffraction lengths a divergent parabolic type diffraction is observed. In nonlinear regime a possibility 
appear: the divergent parabolic type diffraction for large-band pulses to be compensated by the converged parabolic type 
nonlinear Kerr focusing. 
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III. SELF-FOCUSING OF NARROW BAND FEMTOSECOND PULSES. CONICAL EMISSION AND 

SPECTRAL BROADENING 

The laser pulses in a media acquire additional carrier -to envelope phase (CEP), connected with the group-phase 
velocity difference. In air the dispersion is a second order phase effect with respect to the CEP. In linear regime the 
envelope equations contain Galilean invariance, and thus CEP does not influence the pulse evolution. Taking into 
account the CEP in the expression for the nonlinear polarization of third order, a new frequency conversion in THz 
and GHz region takes place. In Laboratory frame, the nonlinear polarization of third order for a laser beam or optical 
pulse, without considering CEP, can be written as follows: 

n 2 E 3 (x, y, z,t)x — xn 2 exp [i(k (z - v ph t)] x 

(18) 

i ^\A\ 2 A + i exp [2i(k (z - v ph t)]A 3 I -I- xc.c, 



while in Galilean coordinates (z! = z — v gr t;t' = t) the CEP, being an absolute phase [2l|, is present in the phase of 
the Third Harmonic (TH) term 

n 2 E 3 [x, y, z,t)x = xn 2 exp [i (k (z' - (v ph - v gr )t')] x 

(19) 

|^|yl| 2 A+ iexp[2i {k (z' - (v ph - v gr )t')]A 3 ^j + xc.c. 

Note that we transform the TH term to a frequency shift of Lu n i = 3ko(v p h — v gr ) = 93GHz in air of the carrying 
wave number Ao = 800nm. The nonlinear amplitude equations for power near the critical one for self- focusing in 
Laboratory and Galilean frame are: 
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respectively. We use AE equations (|20|) and ([21]) to simulate the propagation of a fs pulse, typical for laboratory-scale 
experiments: initial power P = 2Pk r , center wavelength A = 800 nm, initial time duration to = 400 fs, corresponding 
to spatial pulse duration zq = v gr to = 120 /im, and waist rg = 120 y,m. 

Fig. 4 presents the evolution of the spot |A(x,y| 2 of the initial Gaussian laser pulse at distances z = 0, z = 
l/2zdif /, z = Zdiff,z = 3/2zdiff- As a result, we obtain the typical self-focal zone (core) with colored ring around, 
observed in several experiments [Hl-Qjj]. The 3D + 1 nonlinear AE equation (|2"Tj) gives an additional possibility for 
investigating the evolution of the side projection of the intensity \A(x, z'\ 2 profile. The side projection \A(x, z'\ 2 of the 
same pulse is presented in Fig. 5. The initial Gaussian pulse begins to self-compress at about one diffraction length 
and it is split in a sequence of several maxima with decreasing amplitude. Fig. 6 presents the evolution of the Fourier 
spectrum of the side projection k z '\ 2 . At one diffraction length the pulse enlarges asymmetrically towards the 



7 



al a2 




Z=0 Z = 2l^ 

FIG. 4: Nonlinear evolution of the waist (intensity) projection \A(x, y)\ 2 of a 400 /s initial Gaussian pulse (al) at A = 800 nm, 
with spot ro = 120 \im, and longitudinal spatial pulse duration Zo = v gr to — 120 /im at a distance z = Izuff (a2), obtained 
by numerical simulation of the 3D+1 nonlinear AE equation (|20[) . The power is above the critical for self- focusing P = 2P^ r . 
Typical self-focal zone (core) surrounded by Newton's ring is obtained, (b) Comparison with the experimental result presented 
in til- 



(b) 




FIG. 5: fa) Experimental result of pulse self compression and spliting of the initial pulse to a sequence of several decreasing 
maxima [3Q]. (b) Numerical simulation of the evolution of (x,t=z) projection \A(x, z')\ 2 of the same pulse of Fig. 4 at distances 
z — 0, z — Zdiff, governed by the (3D+1) nonlinear AE equation (120[) and the ionization-free model. 




(a) (b) 

FIG. 6: (a) Fourier spectrum of the same side (a;, z) projection of the intensity ^(fca,, k z )\ 2 as in Fig. 5. At one diffraction length 
the pulse enlarges asymmetrically forwards the short wavelengths (high k z wave- numbers), (b) a spectral form observed also in 
the experiments 
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short wavelengths (high wave- numbers). It is important to point here, that similar numerical results for narrow band 
pulses are obtained when only the self-action term in AE equation (f2~Tj) is taken into account. The TH or THz term 
(the second nonlinear term in the brackets) practically does not influence the intensity picture during propagation. In 
conclusion of this paragraph, we should point out that our non-paraxial ionization-free model (|20p and (|21l) is in good 
agreement with the experiments on spatial and spectral transformations of a fs pulse in a regime near the critical 
P > P cr . Such transformation of the shape and spectrum of the fs pulse is typical in the near zone, up to several 
diffraction lengths, where the conditions for narrow-band pulse are satisfied Ak z << ko- 



IV. CARRIER-TO-ENVELOPE PHASE AND NONLINEAR POLARIZATION. DRIFT FROM THZ TO 

GHZ GENERATION 



In nonlinear regime the spectrum of the amplitude function becomes large due to different nonlinear mechanisms. 
The Fourier expression A [k x , h y , ko — k z ,ujo — u>] is a function of arbitrary Ak z = ko — k z and Auj = ujq — uj, which 
are related to the group velocity Aw/ Ak z — v gr (here, we do not include the nonlinear addition to the group velocity 
- it is too small for power near the critical one). Let Ak z denote an arbitrary initial spectral width of the pulse. In 
the nonlinear regime Ak z (z) enlarges considerably and approaches values Ak z (z) ~ ko. To see the difference between 
the evolution of narrow-band Ak z << ko and broadband Ak z ~ ko pulses, it is convenient to rewrite the amplitude 
function in Laboratory coordinates (the dispersion number /3 ~ 2.1 x 10 , being smaller than the diffraction in air, 
is neglected): 

A (x, y, z, t) = B B (x, y, z, t) exp (-i(Ak z (z - v gr t)), (22) 
while in Galilean coordinates it is equal to: 

V (x, y, z', t') = B G (x, y, z' , t') exp {-iAk z z'). (23) 
The Nonlinear Diffraction Equation (NDE) (|20p in Laboratory frame becomes: 



, s ,'dB 1 dB\ A „ 1 d 2 B 
2i(k - Ak z ) — H — I = AB 



dz v gr dt J Vg r dt 2 



n 2 k 2 B 2 1 j \B\ 2 B + i exp [2i ((fc - Ak z )z ~ (k v ph - v gr Ak z )t)]B 3 I + c.c, 



(24) 



and in Galilean frame, the equation (|2Tj) is: 



2(k - Ak z )dG n 1 (d 2 G d 2 G 
= A i G 7T- „ .o — 2.v. 



v gr dt' v 2 gr \dt' 2 " 9r dt'dz' 



(25) 



«2fco B o i 7 |G| 2 G+-exp[2 J ((A :o - Ak z )z' - k (v ph - v gr )t')]G 3 } + c.c, 



where Ak z can get arbitrary values. It can be seen that the nonlinear phases in both coordinate systems are equal 
after the transformation z' = z — v gr t; t' = t: 

{ko - Ak z )z - (k v ph - Ak z v gr )t = (k - Ak z )z' - k {v ph ~ v gr )t'. (26) 

On the other hand, the corresponding frequency conversions are different. In Laboratory frame the frequency conver- 
sion depends on the spectral width Ak z : 

^ni b - k v ph ~ Ak z v gr , (27) 
while in Galilean frame the nonlinear frequency conversion is fixed to the offset frequency 

uj^ 1 = ko{v ph - Vgr) = 31GHz; (A = 800nm) (28) 
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in air. The expression of the nonlinear frequency shift in Laboratory frame (|27[) explains the different frequency arising 
from pulses with different initial spectral width. When the laser is in ns or ps regime, Afc z << fco and the nonlinear 
frequency shift is equal to the third harmonic 3w^ ; ab = 3ojq. In this case (spectral width of the pulse much smaller 
than the spectral distance to the third harmonic), the phase matching conditions can not be met. Thus, the nonlinear 
polarization is transformed into a self-action term. The fs pulses on the other hand have initial spectral width of the 
order Auj? 3 ~ 10 13-14 Hz and for such pulses at short distances in nonlinear regime the condition Aujf 3 ~ <^ ab can 
be satisfied. Thus, the nonlinear frequency shift lies within the spectral width of a fs pulse, and from (|27[) follows 
the condition for THz and not for TH generation. The self-action enlarges the spectrum up to values Afc z ~ fco and 
thus, following (|27|) . the nonlinear frequency conversion in far field zone drifts from THz to ~ 93 GHz [l8j]. Note that 
we consider a single pulse propagation, while the laser system generates a sequences of fs pulses. The different pulses 
have different nonlinear spectral widths when moving from the source to the far field zone. One would detect in an 
experiment a mix of frequencies from THz up to GHz. 



V. NONLINEAR SUB-CYCLE REGIME FOR Ak z « k 



The separation of the nonlinear polarization to self-action and TH, THz or GHz generated terms is appropriate for fs 
pulses up to several cycles under envelope. For fs narrow-band pulses, as mentioned in the previous section, the pulse 
shape is changed by the self-action term, while the CEP frequency depending at the spectral width of the pulse Ak z 
leads to different type of frequency conversion drifts from THz to GHz region. However, when supper-broad spectrum 
occurs (Afc z « fco), the time width of the pulse At becomes smaller than the period of the nonlinear oscillation w^ ; ab . 
In this nonlinear sub-cycle regime, the nonlinear term starts to oscillate with and separation of the self-action 
and the frequency conversion terms becomes mathematically incorrect, due to the mixing of frequencies [321. l33j. For 
the first time such possibility was discussed in 32] , where a correct expression of the nonlinear polarization, including 
Raman response is presented. In the sub-cycle regime the nonlinear polarization at a fixed frequency and Laboratory 
frame becomes: 

n 2 E 3 (x, y, z, t) = n 2 exp [i ((fc - Ak z )z - (k v ph - Ak z v gr )t)]\ x 

(29) 

exp [2i ((fc - Afc z )z - (k a v ph - Ak z v gr )t)]B 3 j , 

and in Galilean frame it is 



n 2 E 3 (x, y, z', t') = n 2 exp [i ((fc - Ak z )z' - k (v ph - v gr )t')]] x 

(30) 



| exp [2i ((fc - Afc z )z' - k (v ph - v gr )t')]B 3 

In spite of the super-broad spectrum, the dispersion parameter in the transparency region from 400 nm up to 800 nm 
continues to be small, in the range of /3 ~ 10~ 4 — 1Q — 5 . The nonlinear amplitude equations for pulses with super-broad 
spectrum in Laboratory system become: 



OA 1 dA\ . . 1 d 2 A 

+ 

(31) 



2*(fco-Afc z) (- + --j=AA--^ + 



and in Galilean frame 



n 2 fcoexp [2i ((fc - Ak z )z - (k v ph - Ak z v gr )t)]A 3 , 



.(ko-Ak z )dv i (d 2 v d 2 v 



v gr dt> x v 2 gr \dt' 2 ar dt'dz' 

(32) 

n 2 fco exp [2i ((fc - Ak z )z' - k (v ph - v gr )t')]V 3 . 
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FIG. 7: Numerical simulations for an initial Gaussian pulse with super-broad spectrum Ak z « fco governed by the nonlinear 
equation (|32[1 , The power is slightly above the critical P = 2Pk r - The side projection |V(a;, z')\ 2 of the intensity is plotted. 
Instead of splitting into a series of several maxima, the pulse transforms its shape into a Lorentzian of the kind V(x,y,z') 
l/[l+x 2 +y 2 + (z' + ia)~ 
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FIG. 8: The evolution of the spectrum \V(k x , k z )\ 2 of the same side intensity projection [T/(o;,a')| 2 
towards small k z wave-numbers (long wavelengths) - typical for Lorentzian profiles. 



The spectrum enlarges 



Fig. 7 shows a typical numerical solution of the nonparaxial nonlinear equation (j3"Tj) (or (|32l) ) for an initial Gaussian 
pulse with super-broad spectrum Afc 2 w fco- It is obtained by using the split step method (4 step Runge-Kutta method 
for the nonlinear part). These results are the same both in Laboratory and Galilean coordinate frames differing only 
by a translation. The side projection |T^(x,z'| 2 of the intensity profile is plotted for different propagation distances. 
Instead of splitting into to a series of several maxima, the pulse transforms its shape in a Lorentzian type form of 
the kind V(x, y, z) ~ 1/[1 + x 2 + y 2 + (z' + ia) 2 + a 2 }. Here, the number a accounts for compression in z' direction 
and a spatial angular distribution. Fig. 8 presents the evolution of the spectrum | V(k x , k z > \ 2 of the side intensity 
projection for the same pulse. The spectrum enlarges forwards the small k z wave-numbers (long wavelengths) - typical 
for Lorentzian type profiles. To compare with Fig. 8, Fig. 9 gives a plot of the side projection \V(k x , k' z \ 2 of the 
spectrum of a Lorentzian profile V(x,y,z') = 1/[1 + x 2 + y 2 + (z' + ia) 2 + a 2 ], a = 2 increases toward the small 
wave-numbers. The numerical experiments lead to the conclusion that a possible shape of the stable 3D + 1 soliton 
can be in the form of a Lorentzian profile. Thus, if we take as an initial condition Lorentzian, instead Gaussian one, 
a relative stability in the shape and spectrum can be expected. Fig. 10 shows the evolution of the \V(x, z'\ 2 profile of 
a pulse with initial Lorentzian shape V(x, y, z', t = 0) = 1/[1 + x 2 + y 2 + (z' + ia) 2 + a 2 ], a = 2. The pulse propagates 
at distance of one diffraction length, preserving its initial shape. 



VI. SPECTRALLY ASYMMETRIC 3D+1 SOLITON SOLUTION 



The numerical simulations in the previous section for broad band spectrum pulses demonstrate a stable soliton 
propagation with a specific initial Lorentzian shape. To find an exact soliton solution, we require that Ak z = kg and 
Alo = wo be satisfied. In air /3 = and the amplitude equation (|3 1 1) can be rewritten as: 
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FIG. 9: Plot of the side projection | V(fca;, k z \ 2 of the spectrum of a Lorentzian profile V(x,y, z) — 1/[1 + x 2 + y 2 + (z + ia) 2 + a 2 ], 
a — 2 increasing towards the small k z wave-numbers (compare with Fig. 8). 




FIG. 10: Evolution of the |V(:r,z'| 2 profile of a pulse with super-broad spectrum Ak z w fco and initial Lorentzian shape 
V(x, y, z' , t = 0) = 1/([1 + x 2 + y 2 + [z 1 + ia) 2 + a 2 ], a — 2, governed by the nonlinear equation (|32[) . The pulse propagates 
over one diffraction length with relatively stable form. 



AB - -jr^ + kln 2 Bl exp [i (2Auj n it)]B 3 = 0. (33) 
v gr at 

To minimize the influence of the GHz oscillation uj n i, we use an amplitude function with a phase opposite to CEP: 

B(x, y, x, t) = C(x, y, z, t) exp(-iAw„;i). (34) 
This corresponds to an oscillation of our soliton solution with frequency u> n i ~ 31 GHz. The equation p3p becomes: 

AC 40 + kln 2 B^ = 2i^^ -*&C (35) 
v z gr at 2 v£ r at v 2 r 

To estimate the influence of the different terms on the propagation dynamics we rewrite equation (|35l) in dimensionless 
form. Substituting: 

t = t Q t; z = z z; x = r x; y = r y; (36) 
r /z = 5~1; z Q = v gr t ; t = 2 x 10" 13 - l(T 14 sec, (37) 

we obtain the following normalized equation: 



(38) 
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where 7 = r 2 k^n 2 BQ is the nonlinear constant, a — 2Auj n ir 2 /v 2 r to and /3 = Au} 2 d r 2 jv 2 gr . For typical fs laser pulse at 
carrier wavenumbcr 800 nm with spot r — 100 /im, the constants of both terms in the r.h.s of equation (|38[) are very 
small (a ~ 10 -2 and /3 ~ 10 -4 ) and can be neglected. Thus, equation (|38]l becomes: 



AC -^+7C 3 = 0. (39) 
Furthermore, we shall assume that the new envelope wave equation Q39p has solutions in the form: 



C(x,y,z,t)=C(P), (40) 

where f — \J x 2 + y 2 + (z + ia) 2 — (t + ia) 2 . From the nonlinear wave equation (|39l) , using (|40[) , the following ordinary 
nonlinear equation is obtained: 



3 dC + d 2 C + 3 _ Q 

The number a counts for the longitudinal compression and the phase modulation of the pulse. When the nonlinear 
coefficient is slightly above the critical and reaches the value 7 = 2, equation (|4"Tj) has exact particle-like solution of 
the form: 



c _ sech(ln(f)) 
f 

Using the fact that exp(ln(r)) — f and exp(—{lnr)) = 4, the solution (|42|) is simplified to the following algebraic 
soliton: 



1 + x 2 + y 2 + (z + ia) 2 - (t + ia) 2 

The solution (|43p gives the time evolution of our Lorentz initial form, investigated in the previous section. As be 
seen from equation (|39p , the solution appears as a balance between the parabolic (not paraxial) wave type diffraction 
of broad band pulse Ak z = kg and the nonlinearity of third order. The maxima of this solution are at the points 
where f 2 = 0. If we turn back to standard, not normalized coordinates, and solve the second order equation z 2 + 
2iaz — 2iav gr t — v 2 r t 2 = 0, only one real solution z = v gr t can be obtained. It corresponds to one-directional 
propagation with position of the maximum on the z - coordinate z — v gr t. As it was pointed above, Fig. 8 presents 
the initial k x ,k z spectrum of the soliton (1431) . While the k Xl k y spectrum is symmetric, the k z projection is fully 
asymmetric, enlarging forwards to low k z wave-numbers (long wavelengths), and has ty pica l Lorentz shape. Recently, 
in experiments with 2 — 3 cycle pulses long range filaments with similar spectral profile [34| are observed. We suppose 
that in this experiment a 3D + 1 Lorentz type soliton was found experimentally for the first time. 



VII. CONCLUSIONS 



In this paper we investigate femtosecond pulse propagation in air, governed by the AE equation, in linear and 
nonlinear regime. The equation allows to solve the problem of propagation of pulses with super-broad spectrum. 
Note that this problem can not be studied in paraxial optics. In linear regime the fundamental solutions of AE (fTJ) 
and DE ^ are obtained and different regimes of diffraction are analyzed. The typical fs pulses up to 50 fs diffract 
by the Fresnel law, in a plane orthogonal to the direction of propagation, while their longitudinal shape is preserved 
in air or is enlarged a little, due to the dispersion. Broad-band pulses (only a few cycles under envelope) at several 
diffraction lengths diffract in a parabolic form. We solve the convolution problem of the diffraction equation DE @ 
for an initial pulse in the form of a Gaussian bullet, and obtain an exact analytical solution (JS|). A new method 
for solving evolution problems of the wave equation is also suggested. We investigate precisely the nonlinear third 
order polarization, including the CEP into account. This additional phase transforms TH term to THz or GHz terms, 
depending on the spectral width of the pulse. Thus, we suggest a new mechanism of THz and GHz generation from fs 
pulses in nonlinear regime. For pulses with power a little above the critical for self-focusing, we investigate two basic 
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cases: pulses with narrow-band spectrum and with broad-band spectrum. The numerical simulation of the evolution 
of narrow-band pulses (standard 100 fs pulses), gives a typical conical emission and a spectral enlargement to the 
short wavelengths. Our study of broad-band pulses leads to the conclusion that their propagation is governed by the 
nonlinear wave equation with third order nonlinear term p9[) . when the THz oscillation is neglected as small term. An 
exact soliton solution of equation (1391) . with 3D + 1 Lorentz shape is also obtained. The soliton appears as a balance 
between parabolic divergent type diffraction and parabolic convergent type of nonlinear self-focusing. Numerically, 
we demonstrate a relative stability of the soliton pulse with respect to the THz oscillations. 
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Fig.l Plot of the waist (intensity's) projection \A(x, y)\ 2 of a 100 fs Gaussian pulse at A = 800 nm, with initial 
spot ro = 60 fiui, and longitudinal spatial pulse duration zq = 30 fim, as solution of the linear equation in local time 
(p~2|) on distances expressed by diffraction lengths. The spot deformation satisfies the Fresnel diffraction law and on 
one diffraction length z = Zdiff the diameter of the spot increases twice, while the maximum of the pulse decreases 
with the same factor. 

Fig. 2 Side (x, r) projection of the intensity \A(x, r)| 2 for the same optical pulse as in Fig. 1. The (x, y) projection 
of the pulse diffracts considerably following the Fresnel law, while the (r) projection on several diffraction lengths 
preserves its initial shape due to the small dispersion. The diffraction - dispersion picture, presented by the side (x, r) 
projection, gives idea of what should happen in the nonlinear regime: the plane wave diffraction with a combination 
of parabolic type nonlinear Kerr focusing always leads to self-focusing for narrow-band (Afc z << fco) pulses. 

Fig. 3 Side (x, z) projection of the intensity \A(x, z')\ 2 for a normalized 10 fs Gaussian initial pulse at A = 800 nm, 
Afc z ~ fco/3, zq = ra/2, and only 3 cycles under the envelope (large-band pulse Afc z ss fco), obtained numerically from 
the AE equation (J5J in Galilean frame. At 3 diffraction lengths a divergent parabolic type diffraction is observed. In 
nonlinear regime a possibility appear: the divergent parabolic type diffraction for large-band pulses to be compensated 
by the converged parabolic type nonlinear Kerr focusing. 

Fig. 4 Nonlinear evolution of the waist (intensity) projection |A(x,y)| 2 of a 400 fs initial Gaussian pulse (al) at 
A = 800 nm, with spot ro = 120 fim, and longitudinal spatial pulse duration zq = v gr to = 120 [im at a distance 
z = 2zdiff (a2), obtained by numerical simulation of the 3D+1 nonlinear AE equation (|20[) . The power is above 
the critical for self-focusing P — 2Pf, r ■ Typical self-focal zone (core) surrounded by Newton's ring is obtained, (b) 
Comparison with the experimental result presented in [12j . 




Fig. 5 (a) Experimental result of pulse self compression and spliting of the initial pulse to a sequence of several 
decreasing maxima [30j . (b) Numerical simulation of the evolution of (x,t=z) projection \A(x, z')\ 2 of the same pulse 
of Fig. 4 at distances z = 0,z = Zdiff, governed by the (3D+1) nonlinear AE equation (l2Tfl) and the ionization-free 
model. 

Fig. 6 (a) Fourier spectrum of the same side (x,z) projection of the intensity ^(fc^, fc z )| 2 as in Fig. 5. At one 
diffraction length the pulse enlarges asymmetrically forwards the short wavelengths (high fc z wave- numbers) , (b) a 
spectral form observed also in the experiments (l2T |. 

Fig. 7 Numerical simulations for an initial Gaussian pulse with super-broad spectrum Afc z ss fco governed by the 
nonlinear equation (f3"2"j). The power is slightly above the critical P = 2P^ r . The side projection \V(x, z')\ 2 of the 
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intensity is plotted. Instead of splitting into a series of several maxima, the pulse transforms its shape into a Lorentzian 
of the kind V(x, y, z') ~ 1/[1 + x 2 + y 2 + [z 1 + ia) 2 + a 2 }. 

Fig. 8 The evolution of the spectrum \V(k x ,k z )\ 2 of the same side intensity projection \V(x, z')\ 2 . The spectrum 
enlarges towards small k z wave- numbers (long wavelengths) - typical for Lorentzian profiles. 

Fig. 9 Plot of the side projection \V(k x , k z \ 2 of the spectrum of a Lorentzian profile V(x, y, z) = 1/[1 + x 2 + y 2 + 
(z + ia) 2 + a 2 ], a = 2 increasing towards the small k z wave-numbers (compare with Fig. 8). 

Fig. 10 Evolution of the \V(x, z'\ 2 profile of a pulse with super-broad spectrum Ak z « kg and initial Lorentzian 
shape V(x, y, z', t = 0) = 1/([1 + x 2 + y 2 + (z' + ia) 2 + a 2 ], a = 2, governed by the nonlinear equation (|52"]h The pulse 
propagates over one diffraction length with relatively stable form. 



